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ABSTRACT 

Analysis of the absolute value of Titan's albedo and its variation with 
increasing phase angle has yielded constraints on the optical properties and 
average particle size of the aerosols responsible for the scattering of visible 
light. The real index of refraction of the scattering material lies within the 
range 1.5 ^ n r ^ 2.0 and the average particle size is somewhere between 
0.2 fim and 0.4 fim. The amount of limb darkening produced by these 
models leads to an occultation radius of " 2700 km. 


The reflection of visible light by Titan is believed to be due chiefly to 
scattering by an optically thick layer of particles analogous to the blue absorbing 
"Axel dust" present in the upper atmospheres of the outer planets. By comparing 
observed properties of Titan with model calculations, we have obtained preliminary 
estimates of some of the properties of the aerosols, including particle size and 
spectral absorption characteristics. In addition, these characteristics allow us to 
predict Titan's limb darkening, thereby removing an ambiguity in the determination 
of Titan's radius found from lunar occultation observations. 

We evaluated the scattering characteristics of our aerosol models with a com- 
puter program based on the doubling method that provides an accurate solution to the 
multiple scattering problem. The single scattering properties of the aerosols were 
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computed using a scheme developed for nonspherical particles (Pollack et al., 1977). 
The free parameters of the model include the real part of the index of refraction n f ; 
the imaginary index n^; optical depth of the aerosol layer at a reference wavelength 
of 0.55 pm, r ; particle size distribution function f(r), where r is the radius of an 
equal volume sphere; and three parameters, a FTB, and SAR, which are related to 
the nonspherical nature of the particles. We used the two parameter size distribution 
proposed by Hansen and Hovenier (1974), oince the scattering properties of the aero- 
sol layer depend almost entirely on one of these parameters; r, the ci'oss section 
(geometric) weighted average particle radius. The second parameter b is a measure 
of the width of the distribution function. The nonspherical parameters a Q , FTB, and 
SAR are, respectively, the ratio of particle circumference to wavelength below which 
the particles act like spherical scatterers and above which they depart from Mie 
scattering; the ratio of light singly scattered into the forward hemisphere to that 
scattered into the backward hemisphere for the non-Mie domain; and the ratio of 
actual surface area to that of an equal volume sphere. 

In the calculations discussed below, we selected the following values for the 
above parameters: n r =1.3, 1.4, 1.5, 1.7, 1.8, 2.0, 2.5, 3.0; n. (A) to be found 
from the observations; r = 10 (essentially infinite optical depth); r to be determined 
from the observations; b = 0.05 (a narrow distribution); a Q = 8 (typical of particles 
lacking sharp edges); FTB = 2; and SAR = 1. 3. 

Even though the phase angle dependence of Titan can only be observed over a 
6.5° range from Earth, it is still possible to obtain some information on the mean 
particle size r from this data. The magnitude of the phase angle variation depends 
on the shape of the single scattering phase function near a scattering angle 9 of 180 . 
This in turn depends on r. When a = 2rrr/A«l, where \ is the wavelength, the aero- 
sol exhibits a Rayleigh scattering phase function, which decreases only slightly with 
decreasing e when 0 is near 180°. When a » a Q , the phase function increases with 

decreasing 9, but when l<a<o: , it decreases significantly as 0 decreases slightly 

c ^ • 
from 180 . Thus, the second a domain exhibits the smallest decrease in brightness 

with increasing phase angle (decreasing 9 ), the third a domain exhibits the largest 

decrease, and the first a domain exhibits an intermediate behavior. Hence, some 

bounds on r can be obtained from the observed phase effect. 
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We have computed the phase effect of our model aerosol laye 'S for a wide 
range of values of r. For each choice of r, we have evaluated m(X) by demanding 
that the computed geometric albedo match the observed value at wavelength X (Nelson 
and riapke, 1978). We compared the predicted phase effects with those observed by 
Noland, et al. (1974) at six wavelengths ranging from 0. 35 jam to 0. 75 pm. The large 
decrease in brightness with increasing phase angle found at the shorter wavelengths 
implies that 1 < a < o- 0 i or 0. 05 pm < r < 0. 45 pm. Within this size domain it is 
possible to obtain an actual value for r from the data. Figure shows a plot of the 
ratio of Titan's disk integrated intensity at 6. 4° and 0° phase angle as a function of 
wavelength. Filled squares and vertical lines indicate the observed values and their 
associated error bars, as found from a least squares fit to the observations, while 
the curves indicate the predicted behavior of several models that come closest to 
fitting the data. The middle curve corresponds to a model with =1.5, r = 0. 35pm, 
while the top and bottom curves refer to r = 0. 32 pm and r = 0.40 pm, respectively. 
We see that the middle curve is not only capable of matching the observed phase 
effect at two wavelengths, as might be expected since there are two free parameters 
n r and r, but is also able to reproduce the observed spectral dependence over all six 
wavelengths. We also see that r can be altered by at most a few hundredths of a 
micron from its optimum value before the fit to the observations becomes unaccept- 
able. However, the value of r can be changed by a somewhat larger amount if n r is 
also varied. 

In Figure 2, the model curves coming closest to fitting the observations are 
shown for several different values of n^. In addition to the cases shown, fits to the 
observations were attempted using n r = 1. 3, 1.4, 2.5, and 3.0. A refractive index 
of 1.4 or less is unable to produce a phase variation much less than 0.98, which 
makes it clearly unacceptable at the shorter wavelengths. Refractive indices greater 
than 2 yield insufficient variation of the phase function with wavelength, a trend which 
can already be detected in the model curve shown for n f = 2.0. So judging from the 
results shown in Figure 2 , 1. 5< n < 2. 0 and 0. 2 pm < r < 0. 4 pm. 

The inferred value of r is consistent with the requirements of the inversion 
model of Titan (Danielson et al. , 1973). This model requires that the mean size of 
the aerosols absorbing some of the incident sunlight be small enough so that they are 
poor radiators in the thermal region of the spectrum, i.e. or « 1 for X ~ 10 pm or 
longer. This condition allows the upper atmosphere to assume higher temperatures 
than the effective temperature. 
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SAMPLE FITS TO PHASE VARIATION 



X, 


Figure I Illustration o] the phase effect calculate J for Titan a ' a fu’htion of uacehngth, for n r — 
1.5 and three different mean particle size i. The ulaes denied from the ohien at ions of Noland cf al. 
( !9~’a i are mdhated by filled squares u ith error />a> i , 
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Figure 3 displays the imaginary index of refraction as a function of wave- 
length, as found by matching Titan's geometric albedo for an assumed visible "sur- 
face" radius of 2700 km. The four curves of this figure correspond to the four models 
shown in Figure 2. We see that the absorption coefficient of the aerosols deci’eases 
by about one order of magnitude between 0. 35 pm and about 0. 6pm, but that it flattens 
toward longer wavelengths. Both the deduced spectral shape and approximate absolute 
value of m should provide useful constraints on the composition of the aerosols. 

Finally, we consider the limb darkening behavior of our most successful 
models. In Figure 4, I(p)/l(p = 1) is plotted as a function of p at wavelengths of 
0.49 pm and 0. 62 pm, where p is the cosine of the angle between the local vertical 
and the line of sight. These predictions refer to zero degrees phase angle. Also 
given in the figure for comparison is the behavior of a Lambert surface. 

These results have relevance for the inference of Titan’s radius from lunar 
occultation observations (Elliot et al. , 1975). These data do not permit the simul- 
taneous solution of both the limb darkening law and the satellite's radius at the 
occultation level. For a uniformly bright satellite (I(p) = I (P = 1)), the oc , . ■>. n radius 
is about 2500 km, while for a Lambert law, it has a value of 2900 km. t 1 l - ed 

limb darkening illustrated in Figure 4 implies that the occultation radius <un oout 

2600 to 2700 km. We plan to make a more precise determination of this important 
quantity. 
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IMAGINARY INDICES OF REFRACTION 



A, /JlTI 


Figure }. Imaginary index of refraction as a function of u^velenglh for the cases shown m Figure 2 
Sumhering of the curies is the Sane as in Figure 2. 



TITAN LIMB DARKENING 



Ftgnre 4- Limb darkening a! uavtleugtht .nmtfmuling ;♦ tun of the chaunelt «W n tht Umar 
acadtalem observations of Filial etll. ( /9'5 J. The numbering of the cnmt if the tame at em Figure 2 
Tht lima darkening of a Lambert surface it shorn far imuparifm. 


DISCUSSION 


D. MORRISON: It seems clear that the surface properties as they could be 
derived from photometric or poiarimetric measurements are important for dis- 
tinguishing between the models. You could have a surface, John, where Don has the 
top of a cloud level. I would like to read to you from Veverka's paper given at the 
Titan Workshop (p. 54) — "the single most important conclusion to be drawn from the 
photometry and polarimetry of Titan is that a Saturn-like cloud model may be required 
to explain the sum of the observations. " Can you comment on this ? 

D. HUNTEN: That is a perfectly reasonable interpretation, but it's not unique. 
On p. 57 of that book 1 point out that Titan could also be paved with a glassy layer of 
asphalt, which v ->uld have polarization properties indistinguishable from those of a 
cloud. And that's not an entirely improbable kind of appearance for a Titan surface. 
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J. CALDWELL; 1 don't think it's quite right to conclude that the cloud has to 
be like Saturr. The optical depth of the atmosphere and dust layer, even in these low 
surface pressure models, is very large in the visible and ultraviolet. It seems to me 
perfectly reasonable that the polarization might be due to the properties of the dust. 

B. SMITH; What optical depth do you use? 

J. CALDWELL: Podolak and Danielson calculate an optical depth due to dust 
of ~ 5 at 5000 X. 

B. SMITH: Either way. there is no hope of seeing the surface. 

D. HI NTEN: And no hope of seeing the cloud top of my model either. 

J. POLLACK: The successful model by Rages and myself shows, starting 
about 0.6 pm, a tremendously sir; rp increase in the imaginary index tc shorter wave- 
lengths; at somewhat longer wavelengths, it tends * flatten out. So, in effect, what- 
ever is making up this layer is something that ha y strong and very sharp 

ultraviolet absorption band. I think that ixjth the shape and the absolute absorption 
coefficient here may allow us to choose between different compositional possibilities. 

With a model we can, of course, define what the limb darkening of Titan would 
be like, and the relevance of that is that the lunar occultation measurement of the 
size of Titan is very dependent upon what limb darkening law you assume. If you 
assume a uniformly bright surface, which would be a flat line in Figure 4, then you 
get a radius of 2500 km. If you assume a Lambert law, you get about 2900 km. 

What our models would say is that you're somewhere in between, but closer 
to the laimbert surface, so therefore, the occultation radius would be something on 
the order of 2600 to 2700 km. 

Finally, 1 would like to note that I’ve become increasingly impressed that the 
apparent secular brightening of Titan is, in fact, a real phenomenon, and there’s now 
something like six years of observations that indicate that Titan’s brightness has 
increased by maybe five or so percent. If we assume that that's right, and we also 
assume that what we sec is a photochemical layer, then it's quite conceivable that 
plausible solar variability could indeed have a very interesting climatic feedback. 

For example, if you have ultraviolet solar variations, which we know from satellite 
observation do occur, variations in intensity will cause a variation in the production 
rate of the smog, which in turn could affect the size of the particles somewhat, and 
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therefore, affect the overall brightness oi Titan, if this sort of linkage is true, it’s 
very interesting in the sense that a very small energy change in the Sun is able to 
enormously amplify the amount of solar energy that is deposited in Titan. 

D. Hl'NTEN: But we should also remember that ordinary seasonal variations 
may be the answer. 

J. CALDWELL: Badolak and Danielson have done a lot of work to derive the 
photometric properties of the dust, and they have a size of about 0. 1 pm. 

J. I\)LLACK: This size (0. 1 pm) is incompatible with the phase angle varia- 
tions (see Figure 1 of this article). Also, they make some arbitrary assumptions about 
the analytic dependence of the imaginary index; they assumed a powc* law. I think you 
can see from Figure 3 of this article that this is not true over live whole spectral range 
of the obs. rvations. 

I>. SMITH: What do you mean by the particle size'.* 

J. POLLACK: Anytime you derive particle size information from brightness 
measurements, what you’re really doing is just deriving one gross property of the 
size distribution, namely, the cross-sectional averaged particle size. For Podolsk 
and Danielson, the number given is the maximum radius of a flat distribution, and the 
effective size is a bit smaller. 

J. POLLACK- On the difference between the equivalent widths on Titan and 
Saturn for the methane bands, it could equally well bo the result of the differences in the 
properties of the scattering medium that's present in the two atmospheres. 

L. TR A 1-TON: IVssibly. But look at the differences between the s| vet rum of 
Jupiter and Saturn where you have appreciably different haze layers. The haze layer 
of ammonia is really thin on Jupiter, it's really thick on Saturn, and yet these big 
differences in the aerosol structure between the two planets lead to only small differ- 
ences in the shapes of the spectra when you compare them against the shapes of 
methane in Titan's atmosphere. 

J. IYH.LACK; Let me be a little bit more specific. In the ease of Jupiter 
and Saturn, there is an optically thin haz.e layer and then a fairly thick, dense cloud 
layer beneath that. In the ease of Titan, there is an optically thick haze layer which 
means that a lot more multiple scattering happens in Titan's atmosphere, while on 
Jupiter and Saturn it’s closer to simple reflecting. 
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1.. TRAFTON: Yes, but there are really three different regimes, 1 really 
don’t see Jupiter anti Saturn being in one regime :md Titan being entirely different. 

But in terms of differences it; the spectrum, l do set' very little differences between 
the shapes of the spectrum of Jupiter and Saturn. 

J. POLLACK; That is what l would expect for the photochemical base; it 
is much thinner in the case of Juuitcr ami Saturn than it is on Titan. 1 believe most 
of the line formation takes place within the actual dust in the case of Titan. 

D. MORRISON: Why should Titan have so much more of this dust or smog 
than Jupiter or Saturn ? 

J. lY'l.l ACK: 1 tnink it goes back to abundances; for one thing, the fractional 
abundance of methane is a Jot more on Titan. 

D. Hl'NTEN: t don’t agree that it's a lot more; what is different is the absence 
of hydrogen in the photochemical processes, so that the reducing power of the atmos- 
phere is negligible. 

J. I\>l LACK: Yes. and nit nogen in Titan’s atmosphere could also be a 
critical clement there. 

1.. TRAIT ON: What is needed is Liberal rv data for methane at a very low 
pressure tike a hundredth atmosphere to maybe a tenth to know whether or not there 
is indeed any pressure dependence at the methane level. 

J. OAl.DWFl 1-: Such a measurement requires a very long path length. 

Lutz aiylied tremendous lengths alremly to do this. You may be asking the 
impossible. 


SUMMARY OF DISCI SSION AMONU TRAFTON, HI NTKN, AND POLLACK: Fink, 
Renner, and Dick (1977> find Im'Iow SOOO A that the bands obey Peer's law and arc pressure- 
independent. Their explanation is that there are many lines per half-w idth even at 
zero pressure. Pressure effects begin to l>o evident at longer wavelengths and might 
be expected in the wings of all bands. lMdolak and diver have explored possible 
saturation effects for the lowest pressures. Temperature effects must also be kept 
in mind. 

M. KLEIN: (addressed to Hunten) If you really have this- cloud at the 77 K level 
of GOO g cm - ”, the opacity to microwaves at 3 mm is no longer ui.appreeiabie. This 
effect may not be insignificant. 
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D. HUNTER: The particles would have to be very large to have a significant 
scattering opacity. Pure methane, either liquid or solid, should have very low 
millimeter-wave absorption, because the molecule is nonpolar. 1 think it would need 
polar impurities to be much of an absorber. 

In addition, 1 don't really believe the cloud is that dense. 1 believe that most 
of that mass must precipitate out and leave a much thinner eioud. All 1 did in the 
paper was follow the Lewis-type prescription in which you condense out all the mass 
that is available at each height and call it a cloud. 

D. STROBEL: Do your methane cloud properties satisfy the observations 
that Low and Rieke tAsirophys. J. 190, 1.143, 1974) made at 5 pm ? It could be 
thermal emission at 103 K or solar reflection with an albedo of 0. 10. 

D. HUNTER: My methane cloud top is about 79 or 80 K, but I don't see why 
it couldn't explain the observations as reflected sunlight. 

J. CALDWELL; Is your cloud opaque at all in the 10 pm region ? There are 
many transparent gaps between the fundamental bands of CM , C H . and C..H , in 
this region. 

D. HUNTER: 1 would think so. In one of my models I did simply postulate 
an opaque cloud. 1 rationalize this by suggesting that Axel dust dissolves in methane 
and gives it some additional absorption that pure methane wouldn't have. 

J. CALDWELL; In my models, the optical depth of the Axel dust is the order 
of 0.05. If the optical depth get- much larger than that, the emission cf the dust will 
go up at 10 pm in contradiction to the observations. 

D. HURTEN; But I'm not talking of dust itself, rather ihe same material afte: 
it’s dissolved in the cold methane clouds. 

J. CALDWELL: My reason for emphasizing this point is ihat any model must 
be consistent with Gillett’s observations at 10 pm, which exclude a high brightness 
temperature. If there's a 300 K surface below the haze, something has to hide that 
from an outside observer. 

D. HUNTER: I think pressure iniueod hydrogen will do that. Even at 10 pm, 
with twenty atmospheres of nitrogen and several km-A of hydrogen at the surface, 
that's a very opaque medium. It's very much, like Venus, 

G. SISCOE: What is the hydrogen escape rate from Titan? 

9 

D. HUNTER: My estimate in Planetar* Satellites is 9 \ It) molecules 
cm ' s . That is based simply on photolysis of the methane we set' to be present, 
and is firm unless someone can find a stronger source. 
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